Nature of perturbation theory in spin glasses 
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The high-order behavior of the perturbation expansion in the cubic rephca field theory of spin 
glasses in the paramagnetic phase has been investigated. The study starts with the zero-dimensional 
version of the replica field theory and this is shown to be equivalent to the problem of finding finite 
size corrections in a modified spherical spin glass near the critical temperature. We find that the high- 
order behavior of the perturbation series is described, to leading order, by coefficients of alternating 
signs (suggesting that the cubic field theory is well-defined) but that there are also subdominant 
terms with a complicated dependence of their sign on the order. Our results are then extended 
to the d-dimensional field theory and in particular used to determine the high-order behavior of 
the terms in the expansion of the critical exponents in a power series in e = 6 — d. We have also 
corrected errors in the existing e expansions at third order. 
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I. INTRODUCTION 

The theory of spin glasses in finite dimensional systems 
has traditionally been approached by the loop expansion 
around Parisi's mean-field replica symmetry breaking so- 
lution 0. However, the picture of spin glasses which 
emerges from this perturbative approach is quite different 
to the one arising in droplet theory 2]. The motivation 
for the present paper was to investigate the possibility 
that perturbation theory in spin glasses might fail for 
some reason (for example, non-perturbative terms like 
"droplets" might dominate their free energy, at least in 
the low-temperature phase). We have started the pro- 
gramme with a study of the nature of the perturba- 
tion expansion in the high-temperature or paramagnetic 
phase, and postpone to another paper the discussion of 
the low-temperature phase. 

In general the nature of perturbation expansions in dis- 
ordered system is far from trivial. This is in contrast to 
the perturbation theory in pure systems, where, for ex- 
ample, the Borel summability of the series leads to an ac- 
curate evaluation of critical exponents . In the study of 
the perturbation expansion for disordered ferromagnets. 
Bray et al. Q found that even in zero dimensions, the 
high-order behavior of the perturbation expansion was 
surprisingly rich. The high-order expansion coefficients 
are sums of two kinds of terms: one type has an unusual 
cosine-like oscillation with increasing periodicity and the 
second type has a simple alternation in sign which domi- 
nates for small disorder. This unusual behavior has been 
further studied in Refs. 0, 0, with the final conclu- 
sion that the series is still summable, but that the simple 
Borel procedure needs to be modified to deal successfully 
with long series. We shall find that the perturbation ex- 
pansion for spin glasses has remarkably similar features 
to those of the disordered ferromagnet. 

Our investigation of the nature of the perturbation ex- 
pansions in the high-temperature phase of spin glasses 



starts by looking at the perturbation expansion of the 
zero-dimensional spin glass problem (which will be re- 
ferred to as the "toy" problem). The zero-dimensional 
field theory is the key to the analysis of the d-dimensional 
field theory as the extension to the d-dimensional field 
theory and critical exponents is a relatively straightfor- 
ward extension of the toy problem . Apart from being 
simple integrals, the zero-dimensional toy field theory has 
the advantage of allowing an analysis without replicas. 
In this paper this is achieved by mapping the problem to 
that of critical finite-size scaling in a modified version of 
the spherical spin glass f3| and this mapping allowed us 
to find the high-order behavior without the use of repli- 
cas. However, to obtain the high-order behavior of the 
perturbation expansion for the d-dimensional field the- 
ory requires the use of replicas and we found that this 
needed the use of a non-trivial replica symmetry break- 
ing scheme in the toy model in order to get results consis- 
tent with our mapping to the spherical model. Our chief 
result is that the perturbation theory is well-defined and 
the dominant high-order terms in the perturbation ex- 
pansion have coefficients of alternating signs. However, 
the perturbation series of the zero-dimensional spin glass 
field theory is not Borel summable in a straightforward 
way due to the presence of subdominant terms. 

The paper is organized as follows. In the next sec- 
tion, we consider the cubic replica field theory of spin 
glasses and obtain the first few expansion coefficients of 
the zero-dimensional toy problem by explicitly evaluating 
the Feynman diagrams. In Sec. IIIII we show the equiva- 
lence of the zero-dimensional field theory and the critical 
finite-size scaling of the modified spherical model. Using 
this mapping we obtain the high order behavior of the 
perturbation expansion for the toy problem in Sec. lIVI In 
Sec.0 we consider the toy problem using replicas. This 
is generalized to the problem of the high order terms in 
the e expansion in Sec. I VII We conclude with a discussion 
in Sec.lVHl 
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II. REPLICA FIELD THEORY OF SPIN 
GLASSES 



The replica field theory of spin glasses, (see Ref. 
for a review), starts from the Hamiltonian density 
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As usual the field components qaf3 {a ^ (3 and qap = qpa) 
take all real values, and the indices such as a take the 
values 1,2, 3,..., n. In the limit when n goes to zero, 
such a field theory is thought to capture the physics of 
finite dimensional spin glasses. The quartic terms work 
as stabilizing terms, but for d < 6 are irrelevant variables 
which we shall drop. One question which we shall ask is 
whether the resulting cubic field theory is well-defined in 
the limit n — > 0. It is possible this approach may not 
be valid as cubic theories have Hamiltonians which are 
not bounded below 0]. However, there are examples of 
field theories existing where unphysical limits, such as 
the number of field components n is taken to zero, which 
saves these apparently unphysical theories 0, E| ■ This 
seems also to be the case for spin glasses as our work 
shows that the coefficients in the perturbation expansion 
in w alternate in sign, which is an indication that the field 
theory remains well-defined in the limit when n goes to 
zero. 

The partition function of the zero-dimensional spin 
glass field theory is given by 



n 

Q</3 



dqap 



exp 
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a,f3 



+ 2^ qapqpjq^a 



(2) 



The perturbation expansion is well-defined irrespective 
of whether the integral of Eq. (O actually exists. The 
perturbation expansion in w yields a series 
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(3) 



where we take — vP' j (t/2)'^ as the expansion parame- 
ter of the problem. The use of t/2 instead of r is for later 
convenience. The series expansion of the corresponding 
free energy is given by 
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(4) 



Although this zcro-dimcnsional theory is nothing but 
a multiple integral, it is not an easy task to calculate the 
expansion coefficients Ak or Bk directly by expanding 
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FIG. 1: Feynman diagrams to 0{g ) contributing to the 
renormalized propagator with their corresponding contribu- 
tions. 



the exponential in ([SJl, because of the complicated struc- 
ture of the internal summations. In fact, even using a 
symbolic manipulation program on a computer, we find 
it very difficult to get the expansion coefficients higher 
order than the first couple of terms. Fortunately, in the 
study of a cubic field theory similar to ours - the perco- 
lation problem ~ different types of the internal con- 
tractions occurring in the theory were classified diagram- 
matically, and all the relevant diagrams were given up to 
0{g^). To this order there are five different internal con- 
traction types. (See the diagrams denoted by a, (3, 7, 5 
and A in Figs. and |21) These were translated into the 
spin glass problem in Ref. In Figs. ^E^ndl^l all the 
diagrams to this order contributing to the renormalized 
propagator and vertex are listed along with their contri- 
butions. Note that, in the zero-dimensional field theory, 
each contribution is just given by the product of these 
contraction factors. 

The expansion coefficient Bk of the free energy can be 
obtained by noting that —d{PF)/dw is just the renor- 
malized three-point function, which can be written in 
terms of the renormalized propagator and vertex. Col- 
lecting the contribution from each diagram and using the 
results of the internal contractions, we calculate the per- 
turbation expansion coefficient Bk up to four-loop order 
(0(58)). We find that Bk = \n{n - l){n - 2)//f/8^, 
where /i = 1/2, /z = n- 2, /g = - 31ri/2 44/3 and 
/4 = 22n^ - 123^2 + 229n - 148. Therefore, the free en- 
ergy for the toy spin glass problem, lim„^o F /n, is given 
by a power series in with coefficients of alternating 
signs up to 0(5*). 

In the course of investigation, we have discovered errors 
in two {5 and A in Figs. H ISJ of the five types of 
contraction reported in Ref. |l2| . The correct results we 



obtain are 5 



9n^ + 54n - 104 and A = bn'' - 14n, 
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FIG. 2: Feynman diagrams to O(g^) contributing to the 
renormalized vertex with their corresponding contributions. 



while the results in Ref. were 6 = — 3n^ + 38n — 94 
and A = 5n^ — 2n — 12. The latter can easily be shown 
to be incompatible with the simple n = 3 case. These 
corrections will change the 0(e'^) terms in the critical 
exponents in d = 6 — e dimensions. From the explicit 
expressions for the critical exponents 77 and v obtained in 
Ref. ^3 1 calculate the correct form of the e expansion 
to third order to be 



= -0.3333e+ 1.2593e2 + 0.7637e3 (5) 
- 2 + ri^-2e + 9.2778e^ - 6.4044e^ (6) 



The corrected series for — 2+77 shows an oscillation in 
signs in contrast to the one in Ref. where the 0{e^) 
term was positive. 



III. MAPPING OF THE TOY PROBLEM ONTO 
A MODIFIED SPHERICAL MODEL 

As mentioned in the Introduction, we study the zero- 
dimensional field theory by mapping it onto a modified 
version of the spherical spin glass model. The Hamilto- 
nian of the spherical model is 



— „ JijSiSj 



(7) 



2J 



with the spherical constraint J^i — ^ among the N 
spins. Unlike the conventional spherical spin glass model 
, we take Ju ^ in addition to the infinite-ranged 



interactions Jy = Jji {i / j). They are chosen from 
Gaussian distributions 



e 2 Jij . 



(8) 



The presence of the diagonal interaction does not make 
any difference in the N ^ 00 limit. In the following, 
however, we consider finite size corrections in this model. 
The partition function can be written as 

^ f /:n-./::s 

[^{Nz-zY,Sf+J2j.,S,S,)], (9) 



X exp 



where (3—1 /T is the inverse temperature and the chemi- 
cal potential z was introduced to represent the delta func- 
tion SiN-j:,Sf). 

In order to make a connection to the replica field the- 
ory |(2J), we replicate the partition function n times and 
average over the Gaussian bond distribution. We then 
take the usual Hubbard-Stratonovich (HS) transforma- 
tions on the factor exp[{/3V47V) ^^.^(E^ S^Slf] to get 
the spins on the same site, where the Greek indices de- 
note the replica components. In order to do that, we 
need to introduce the diagonal qaa and off-diagonal qa-y 
{a ^ 7) fields for the corresponding factors. We have 



exp 



X exp 



(10) 



and 



exp 



£1 
2N 
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a<7 



(11) 
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We can then integrate over the spin variables to obtain 
the replicated partition function as integrals over the HS 
fields, Qaa and Qaj and over the replicated chemical po- 
tential z„: 



N 
47r/32 



a<7 ^ 



n dza r N ^ 2 N sr 2 

■(^/3z„-trln [(/? 

\)l (12) 

a 

where tr is taken with respect to the replica index. 



N 
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In the large- iV limit, these integrals can be evaluated 
by the steepest descent method. For T > = 1, the sad- 
dle points are at q^-f = 0, qaa = 0^ and = 13 + 
This is the well-known result 8] from the N ^ oo 
analysis of the spherical spin glass. We investigate the 
finite-size corrections in the limit T ^ Tc = 1 by con- 
sidering the fluctuations around these saddles. Writing 
Qaa = 0^ -\- Va and I3za = 1 -\- (3"^ + ixa, we have 



N 



^{Va - iXa) 



(13) 



- ytrln [{1 - {ya - iXa)}5a-f - qa-r] 



where C = exp((nA^/2)(l -1-/3^/2)). If we expand the log- 
arithm in powers of the fields, we find that the quadratic 
terms in the diagonal fields Xa and ya inside the expo- 
nential are given by 



^^Y,^{T^ ^l)yl + 2ix^y^+xl 



(14) 



Diagonalizing this quadratic form, we find two eigenval- 
ues with nonvanishing negative real parts at T = Tj., 
which implies that the diagonal fields are hard modes 
near and can be integrated away without encounter- 
ing divergences. Therefore the critical behavior is de- 
scribed by the off-diagonal partition function, which can 
be written as 



ioff 



n (^7^) "'^^aiexp 



C(<7 



N 



2 

a7 



N 



X! (Ia(3qi3-,q-,a + 0{q'^) 



a,/3,7 



(15) 



where the quartic and higher order terms in g^-^ all have 
coefficients proportional to N. A key point of this discus- 
sion is to note that, in the limit where t = [T — Tc)/Tc 
and N oo, the quartic and higher order terms can be 
neglected if we keep Nt^ finite. This can be easily seen 
by rescaling — > qaj/ VNt in H15|) . In fact we can show 
that the off-diagonal partition function in this limit is ex- 
actly the same as the zero-dimensional cubic field theory 
defined in (j^J and (O after identifying the expansion pa- 
rameter as = l/{Nt^) and r = - 1 = t{2 + t) ^2t 
as t ^ 0. We have Zos{N,f3) Z{j^) as TV oo, 
t ^ and Nt^ — > finite. A similar observation was made 
for the critical finite size corrections of the Sherrington- 
Kirkpatrick model in Ref. (In this paper, the series 
for the free energy of the toy model was given to order 



FIG. 3: The integration contour Csp used in Eq. I|16|l . The 
filled circles represent the eigenvalues J\ schematically. The 
largest and the second largest eigenvalues are denoted by Ji 
and J2, respectively. 



IV. HIGH-ORDER BEHAVIOR OF THE 
EXPANSION COEFFICIENTS: TOY PROBLEM 

A. Leading-order behavior 

Having established the equivalence of the modified 
spherical spin glass model to the cubic replica field the- 
ory, we can analyze the toy problem ^ without using 
the replicas. Integrating over the spin variables in lO, 
we obtain 



^sp = /3 / ^ exp 



47ri 



Nf3z 1 



A 



2 



2 2 V 2 

(16) 

where J\ denotes the eigenvalue of the matrix Jij . Note 
that the contour Cgp of integration lies to the right of the 
largest eigenvalue Ji. (SeeFig. 

For large N, the integral is dominated by the saddle 
point determined by 



N ^ z~Jx 



(17) 



In the limit N 00, one can evaluate the sum on the 
right hand side using the Wigner semicircle law for the 
eigenvalue density p{J\) Here we are interested in the 
finite N corrections, in particular, the limit where — > 
00 and t with Nt^ held fixed. In this case, one has 
to find a solution z of (|17|1 for given bond realization Jij , 
then calculate the free energy —TlnZsp from (|16|l . which 
has to be averaged over the bonds. This is an extremely 
difficult task to carry out analytically. Instead here we 
make a self-consistent approximation where we assume 
the saddle point is located well away (in the sense to be 
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FIG. 4: (a) The disorder distribution responsible for the dom- 
inant high-order behavior. The largest eigenvalue Ji is far 
away from the saddle point that the contour passes through. 
The rest of the eigenvalues (the shaded region) are approxi- 
mated as a continuum distribution, (b) The disorder distri- 
bution for the sub-dominant high-order behavior. The sepa- 
ration A between the eigenvalues Ji and J2 is large so that 
the eigenvalues smaller than J2 can be approximated as a 
continuum. 



specified below) from the largest eigenvalue Ji . Since the 
eigenvalues can be regarded as one-dimensional electric 
charges interacting logarithmically ^14. 15j, an analogy 
to the multipole expansion in electrostatics suggests that 
we can treat the eigenvalues less than Ji as a continuous 
distribution given by the semicircle law. That is to take 



piJx) 



2{N -1) 



Jl-Jl+6{Jx-Ji). 



(18) 



In the Umit of TV -> cx), t = (/3~ 
(Tc = Ji/2), and finite r/ = N^/^t, 
tance of the saddle point from Ji 
OiN-"^/^). In terms of C = N'^'^{z 



1 - Ji/2)(2/Ji) ^ 
we find that the dis- 
is scaled as z — Ji ~ 
- Ji)(2/Ji), the sad- 



dle point equation (|17|l reduces to 77 = ~ 1/C i^i this 
approximation. As will be shown below, the high-order 
behavior of the perturbation expansion is described by 
the small (or large 77) behavior of the partition func- 
tion. The distance of the saddle point from the largest 
eigenvalue (measured in terms of C,) will be large for large 
•q and therefore the approximation can be justified for de- 
termining the high-order coefficients of the perturbation 
expansion. (See Fig. 0] (a).) 

Within this approximation, the sum in Ijltill consists of 
the term involving Ji and 



^^'ln(z- Ja) ~ ^ 



1 

2' 



(19) 



where the prime indicates the largest eigenvalue is ex- 
cluded from the sum. By changing the integration vari- 
able to C in 1)16(1 and taking the large- limit, we obtain 



-N(3fo ^ 

e 6 



A3 



dC_ exp(- lr;C+K^) 
47ri 



e 3 



(20) 



where /3/o = -/?Ji/2 + (1/2) ln(/3Ji/4) + 1/4 -t- and 
the contour C starts from |^| — 00 with arg(^) — —tt/A 
and extends to |^| = 00 with arg(^) = 7r/4. We can 
evaluate the above integral explicitly as 



PN-^e-'^Pf'' efeAi(^) 



(21) 



where Ai is the Airy function. We note that the above 
expression is valid for both T > Tc [rj > Q) and 
T < Tc (ry < 0). Since the argument of the Airy func- 
tion is an even function of 77, we obtain, in the large- 
N limit, Ai{r]^/4) ^ exp[-\ri\^ /12]/ y^2TT\rf\ using the 
asymptotic behavior of the Airy function. Therefore, 
above Tj,, the leading contribution to the free energy den- 
sity, —N~^ In Zgp, in the large- A^ limit is just (3fo. On the 
other hand, below Tc, it is given by /3fo — /6. One can 
explicitly check that these quantities coincide with the 
free energy densities given in Ref. for the spherical 
model when T approaches Tc from above and below. 

We now consider the finite size corrections above Tc- 
To do this we introduce Zi by writing (|20|l as 



(3e- 



-Nf3fo 



V2TiNt 



(22) 



Note that the square root in the denominator of 1)22(1 
comes from the Gaussian fluctuations around the large- 
A^ saddle point. The finite size corrections relevant to 
the zero-dimensional cubic field theory comes from Zi, 
since we can show that Zi — Zi{g^) is a function of the 
expansion parameter — r\^^ only. It is given by 



Zi{g') 



du 
27ri 



exp 



u 
~2 



(23) 



where u — g^^^i- Although Zi{g^) can be evaluated ana- 
lytically as in ((21(1 . we calculate the high-order expansion 
coefficients of Zi by an indirect method, where the self- 
consistency of our approximation is more apparent. (Re- 
member that we are striving to get the behavior of the 
high-order coefficients exactly; our approximation of the 
distribution of the eigenvalues as a continuum described 
by the semicircle law does not give the correct low-order 
coefficients in the expansion of the free energy). From 
the integral representation ((23(1 . we can analytically con- 
tinue Zi(g'^) to any complex g^ by rotating the contour 
appropriately. We find that Zi has a branch cut along 
the 5^ < axis. The imaginary part of Zi is discon- 
tinuous crossing this axis. For small jgp, we can eval- 
uate the discontinuity by the steepest descent method. 
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FIG. 5: The rotated contours used in Eq. 1)23^ for (a) 

<^^S,{9^) = ^ s-nd (b) arg((;^) — —n. The insets show how 
these contours must be deformed to pass through the saddle 
points u — and u = 1 (filled circles). 




FIG. 6: The integration contour in the space used to deter- 
mine the expansion coefficients qk in Eq. The integral 
reduces to the one along the branch cut on the < axis 
as the contribution from the circle vanishes when the radius 
gets large. 



For arg((/^) = ±7r, the contour C is rotated by ±7r/3 as 
shown in Fig. |31 Among the two saddle points u ~ and 
u = 1 that C can pass through, the latter produces a 
real quantity which is the same for arg(g^) = ±7r, while 
the former is responsible for the discontinuous imaginary 
part 

ImZi(52;arg(g2) =±vr) =T^exp(^)[l + 0(52)]. 

^ (24) 

This can be used to extract the coefficients Ak of the 
perturbation expansion. We follow the standard proce- 
dure 16, 17] by writing a dispersion relation for Zi{g^) 
for > in terms of an integral over a contour that 
runs around the cut in the negative axis. (See Fig. El) 
Therefore, the coefficients qk = lini„^o Ak/ti of the per- 
turbation expansion in the toy spin glass field theory is 
given by 



ax 



dg 



2 lmZi{g^;avg{g 



(9 



2\K+1 



(25) 



For large K, this integral is dominated by the saddle 
point = — 1/ {6K) . This implies that the information 
on lmZi{g'^) for small g^ can be used to obtain Ak for 
large which justifies the present approximation. We 
finally obtain 



The coefficients with alternating signs are consistent with 
the low-order behavior obtained in Sec. HII 



B. Subdominant behavior 

There can be other contributions to the free energy 
from different disorder distributions. For example, when 
the eigenvalues are distributed in such a way that the 
saddle point is not very far from the largest eigenvalue, 
the above approximation breaks down. In this case, we 
expect a different behavior of the free energy. To handle 
this, we again make an approximation which can be justi- 
fied self-consistently. We consider a disorder distribution 
where the second largest eigenvalue J2 is well separated 
from Ji such that the spectrum below J2 can be described 
by the semicircle law. (See Fig. 0](b).) This assumption 
is justified in the following analysis. We approximate 



2{N -2) 



ttJI 



J? 



i=l,2 



and 



dz exp - f ln(f ) - i ln(^ - ^a)] 



47ri 



(28) 

where the double-primed sum excludes A = 1 and 2. This 
sum can be evaluated as in HI9|I with J2 replacing Ji. 
We make the same series of integration variable changes 
leading to (EOJ and ^ (using C = N'^l^{z - J2)(2/J2) 
and Tc = J2/2 in this case), and take the large- limit. 
We obtain Zsp ~ (ie~^^^° Z2 where /o is the same as 



(26) before with Ji replaced by J2 and 



du exp[4(^-l^-f i)] 



(29) 
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with the eigenvalue spacing A = N'^/'^{Ji — J2)(2/J2). 
(Recall = ff ) 

The contribution from this arrangement of disorder to 
the free energy, which we denote by -Fsubi is obtained 
by averaging — In Z2 over the distribution p( A) of the 
eigenvalue spacing A. Among the subdominant contribu- 
tions to the high-order behavior, we focus on those from 
possible zeros of Z2. By explicitly evaluating the con- 
tour integral l|29|l numerically for given we find that 
there exists a complex conjugate pair of zeros, Aq and 
Aq in the complex-A plane. To make analytic progress 
on the contribution from these zeros, we look at the 
fluctuation around the saddle point u = 1. By writ- 
ing It = 1 + i.gy and neglecting 0{g) terms, we obtain 
Z2 ~ (Vff72)^2(.9',A), where 




Z2{g\^)^Z2{v)= ^ 



dy 



2tt y/v + iy 



(30) 



with u = (l-g4/^A)/(2g). Note that V is assumed to be 
of 0(1), which means A ~ Olg^^^"^). This is consistent 
with the present approximation where the separation of 
eigenvalues A is very large. Writing 



ds 



-s (v+iy) 



+ iy ^J2-K 
and integrating over y in ()30|l . we have 



ds 
2^ 



: — Ki(—). (31) 



The zeros of Zi come from the infinitely many zeros of the 
modified Bessel function iiTi We can arrange them 

as complex conjugate pairs, Vm and w^, m — 0,1,2,..., 
given approximately for large m by 



e^"/2[-iln2+ (4m + 3)7r] 



(32) 



and for m - 
9.4244- 







exactly (up to four decimal places) by 
0.6928i. Among the infinitely many zeros, 
we focus only on the pair vq and Wq closest to the origin 
as our numerical evaluation of zeros of (|29|l is consistent 
with Ao = (1 — 2gvi^lg'^l'^ for small g. The other zeros 
probably give only subdominant contributions compared 
to the first ones. 

The subdominant contribution to the expansion coef- 
ficient denoted by Hk can be calculated as 



FIG. 7: The deformed contour in g space that leads to 
Eq. 13511 . The branch points go and Qq^ correspond to the 
zeros of the logarithm. 



The integral over g^ is done by deforming the contour 
such that it runs along a circle of radius i? and around 
the branch cuts associated with the zeros, g\ and g'^ of 
Z2 for fixed A. (See Fig.[7|) The integral along the circle 
vanishes as i? — s- 00 and we have 



aK 



dA p(A) 



dg2 



9o 



dA p(A) 



1 



dg' 



9o 



2\K+1 



(35) 



The eigenvalue spacing distribution is known to take 
the form p(A) ~ exp(-(7rVl6)(A/D)2)) for large A 
when the mean eigenvalue spacing is D. From the numer- 
ical results in Ref. 18] on the eigenvalue spacing distri- 
bution of real symmetric matrices, we can calculate the 
mean spacing near the edge of the spectrum as D ~ 2.30. 
The integral over A in (|35|l can be done using the steep- 
est descent method for large K. The saddle point for the 
first term in l|35|) is given by 



8D^ go ydAJ ' 



(36) 



where dgo /dA can be obtained from the defining equation 



aK 



1 
27ri 



dg^ 



2\K+1 



sub ; 



(33) 



where the integral is over a closed contour surrounding 
the origin. By interchanging the order of integration, we 
can write 



aK = I dAp{A) 



ig') 



'^^■^InZ^iv). (34) 



1 - 50 ^ = 2.90^0 • 



(37) 



We can solve (|36ll and H37|) for large K to obtain the 
saddle point values 



A 



sad 



1 - 
1 - 



^4q;8 
2(1- 



5 s 
6a' 
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where a = n'^ / {12D'^) ~ 0.155. Note that the large-if 
behavior corresponds to smaU g and large A with A ~ 
g~^/^, which means that these results are entirely in the 
regime of the present approximation. Inserting these into 
(|35|l we finally obtain 



(5 1 
cK^ + 0{K^) 



(38) 



where a = a''^!^ ~ 4.05, h = -3a3/®Re(wo) 3.15 and 
c = 3Q;3/®Im(?;o) ^ 3.38. Compared with J^HIl, this is 
a subdominant contribution containing only a fractional 
power of X!. The coefficients do not alternate in sign 
as in H26|l but oscillate with a cosine function with an 
increasing periodicity. The situation is similar to that 
in the zero-dimensional disordered ferromagnets 
where this type of oscillation also occurs in the subdom- 
inant terms. We expect that as in the case of disordered 
ferromagnets the subdominant terms make the resumma- 
tion of the series non-trivial such that a straightforward 
Borel summation is spoiled. However, it seems likely that 
the series could be resummed in other ways and anyway, 
the evidence from Ref. Q suggests that the straightfor- 
ward Pade-Borel method works well for short series even 
in the presence of subdominant terms. 

There are obviously other types of bond distribution 
which could give rise to subdominant contributions to the 
high-order behavior of the perturbation series besides the 
one studied in this subsection. We suspect that the type 
studied here provides the largest of these contributions 
but we have no proof of this. 



V. REPLICA APPROACH TO HIGH ORDER 
BEHAVIOR IN THE TOY PROBLEM 

While the mapping of the toy integral, Eq. Q, to 
the spherical model has enabled us to obtain direct esti- 
mates for the high-order terms of its perturbative expan- 
sion, in order to obtain high-order estimates for the d- 
dimensional field theory and hence for critical exponents 
we have to discover how to obtain the same high-order es- 
timates directly from the integral in the replica variables 
qa(i- Once this has been done the extension to field the- 
ory is relatively straightforward and is carried out in the 
next Section. Unfortunately the direct replica approach 
is neither obvious nor rigorous. Without the results ob- 
tained from the mapping to the spherical model we would 
have had no confidence in the replica procedure which we 
were forced to use. 

It is useful to first examine the integral for the special 
case of 71 = 3 when the integrals can be done explicitly 
and exist if w is pure imaginary. This case was analyzed 
in Ref. Here we follow the same analysis to clarify 
some points which will be important in the case of general 



n. For n = 3, we have after setting r = 1 for simplicity 
47r 



dRR^e~ — f{wR^), (39) 



where 



f{wR^) 



47r 



2Tr <.7r 



d6 sin 9e 

Jo 

'3 \ /„,.T?3 



wR sin 9 cos 6 cos (h sin < 



3 5 V6V3/ V6V3/ 



(40) 



with the modified Bessel function I^{x). For pure 
imaginary coupling w — iw' , w' real, / = 
{TT/3)Ji/G{w'R^/6V3)J_i/e{w'R^/6V3) with the ordi- 
nary Bessel function J,y{x). We can easily see that the 
integral in H39() is well-defined and real in this case. 

When analytically continued to real coupling w, Z3 de- 
velops an imaginary part which is discontinuous crossing 
the branch cut along the real w axis. For small values 
of coupling, one can calculate the discontinuity using the 
steepest descent method on the integral in Since 
the main contribution to the imaginary part comes from 
the saddle point R ~ 0{l/w), we first study the asymp- 
totic behavior of / when is very large, which is 



f{wR^) 



V3 



w 



i?3 



fwR^\ 



cxp 



wR^\ 



(41) 

The above expansion is valid for < arg(w) < tt, or for 
< arg(w^) < 27r. For arg(w) = 0, only the first ex- 
ponential in (|41|l is important and the integral in (|39|) is 
dominated by the saddle points Ri — and R2 — VS/w. 
The steepest descent direction at R2 is perpendicular to 
the real-i? axis along which the integral generates the 
imaginary part. We can deform the contour in H39|l such 
that it starts from Ri along the positive real-i? axis to- 
ward i?2 and makes an upward turn at i?2. (Examples 
of similar deformations are in Ref. |0|.) We note that, 
since we only pass a half of the steepest descent path 
of R2 in this way, the Gaussian integral coming from 
the fiuctuation around R2 produces a half of the total 
fiuctuation contribution. Keeping this in mind, we eval- 
uate the integral to obtain lmZ3{w) = exp(— 2^) for 
arg(w) = 0. For arg(it;) = tt, the saddle points are R'l = 
and i?2 = —y/S/w and the continuation of the contour 
used for arg(it;) = to this case is the one which makes 
the downward turn at i?2. We finally obtain the discon- 
tinuity in the imaginary part along the branch cut on the 
real-w axis as 



ImZsiw) = ±exp(--^), 
2w^ 



(42) 



where the positive and negative signs correspond to 
arg(?«^) — and arg(w^) = 2tt, respectively. This ex- 
ponentially small imaginary contribution to the partition 
function at the physical coupling can be used to deter- 
mine the large order behavior of the perturbation expan- 
sion, but at the same time its presence indicates that the 
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cubic field theory for n = 3 is ill-defined and requires 
stabilizing quartic terms for its existence. 

We now study the case of general n. For small values of 
the expansion parameter, w^/(r/2)'^, we can evaluate the 
integrals in Eq. ^ using the steepest descent method. 
Saddle points are found by solving 



(43) 



The trivial solution of this equation qap = is the 
starting point of the perturbative expansion. Non- 
perturbative terms arise from its non-trivial solution. 
The set of saddles which we study are Qaf} — q, when 
both a and /3 lie in the interval between 1 and r, where 
and q — T / {r — 2)w. This can be described schematically 
as 



9a/3 — 














r 


n— r 



(44) 



Clearly other types of saddle exist besides this. We first 
focus on this class of saddles, which we call scheme I. The 
Hamiltonian density for this saddle is given by 



"(r — 1)7 



12(r 



(45) 



There are more solutions to Eq. H43|l with the same 
Hamiltonian H45|l . which can be obtained by switching 
the signs of some of qap- One can take, for example, 
qip = —q for (3^1, and q^p — q for a,j3 7^ 1. Since 
one can pick any other subscript than 1, the number of 
such solutions is r. One can also switch the signs of two 
different sets such as qi^ = (727 = —q for 7 7^ 1,2 and 
keep all the other elements qafj — q. The number of such 
solutions is r(r — 1)/2. Similarly, one can switch the signs 
of 3, 4, . . . , r — 1 different sets of qap. However, we can 
show that the solutions obtained by switching the signs 
of k different sets are equivalent to those from picking 
1 — k different sets for fc = 1, 2, . . . , r — 1. Therefore, the 
total number S{r) of the solutions with the Hamiltonian 
(gni is 



[r/2] 



-yr—l 



2[{t/2] 



r odd, 
^, r even. 



(46) 



where [x\ denotes the largest integer less than or equal 
to X. 

We need to include the Gaussian fluctuations around 
these saddles. The matrix of the second derivatives has 
six distinct eigenvalues: a "breather" mode with eigen- 
value —T which is non-degenerate, one other negative 
eigenvalue, —T/{r — 2) which is (n — r)-fold degener- 
ate, and eigenvalues 2T/(r — 2), (r — l)-fold degenerate. 



rr/(r-2), (r(r- 3)/2)-fold degenerate, (2r- 3)T/(r- 2), 
(r— l)(n — r)-fold degenerate, and r, {n — r){n — r—l)/2- 
fold degenerate. Setting t = 1 for simplicity and collect- 
ing the contributions from the saddle points to Gaussian 
order, we obtain 



Zm(u;2) = [l + 0(«;2)] 
1 " 

+ -^"a5(r-)exp 



(47) 



r=3 

r~2 
2 

r-2 
2r-3 



-{r- 1) 



r 



12(r- 2)2, 

(-(^-2)) 



where the first term corresponds to the usual perturba- 
tion expansion in w. The sum over r starts at 3 since in 
Eq. (|43|l the index 7 must differ from both a and (3. The 
factor "Cr denotes the number of ways of introducing r 
non-zero blocks. Collecting the contributions from the 
saddles in this way, instead of deforming a contour in a 
multi-dimensional complex space, determines the parti- 
tion function up to an overall factor. The analysis of the 
n = 3 case suggests that there is an overall factor of 1/2 
coming from the fact that, for the nontrivial saddles, only 
a half of the Gaussian integrals contribute compared to 
the perturbative one. Indeed, for n = 3, one can explic- 
itly check that the non-perturbative part of (|47|l reduces 
to Eq. H42|l . The negative eigenvalues are responsible for 
the factor (— l)("~''+i)/2^ which generates the imaginary 
part in Z . 

For finite n, the saddles in the scheme I correspond 
to the partition function which is well defined except for 
real w resulting in a branch cut on the positive ■uP' axis. 
The discontinuity of the imaginary part across the cut 
is exponentially small ~ cxp(— n(?i — l)/(12(7i — 2)'^'w'^)). 
(The imaginary part of H47|) is dominated by the r = n 
term, since the Hamiltonian H45|) decreases monotonically 
as r increases.) If the cubic spin glass field theory is well- 
defined for real coupling w, we expect that the cut moves 
to the negative axis as we take the n —> limit and 
that there is an exponentially small discontinuity across 
the negative w'^ axis. The migration of the cut as the 
analytic continuation of rt ^ is taken is exactly what 
happens in the percolation problem 0. In that case, the 
Hamiltonian Tis for the saddles depends explicitly on n 
such that it changes sign as n ^ 0. In the present case, 
Tls in H45() is independent of n, and there is no way of 
producing an exponentially small discontinuity across the 
negative axis from these saddles. Therefore, we con- 
clude that the saddles in the scheme I are not sufficient 
to describe the partition function in the n — > limit. 

This observation leads us to consider another type of 
solution of Eq. (|43|l . which we call scheme II. It is inspired 
by the replica symmetry breaking scheme used by two of 
us to describe the free energy fluctuations. Instead 
of taking all nonzero for a, (3 = l,2...r as in the 
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previous scheme, we set gap = q only when a and jS 
belong to the m blocks of size r/m on the diagonal of the 
r X r matrix as 



\ 



(48) 



The key point of the construction of these saddles is 
that we let m — > cxD before rt — > 0. For finite m, this 
scheme is just a generalization of the scheme I which 
produces only subleading terms. For m — > oo, however, 
the Hamiltonian becomes 



Hs = 



48^2' 



(49) 



which has the opposite sign to the one in the scheme I. 
This can now describe the partition function where the 
cut lies on the imaginary w axis. The solution to H43|l is 
q = T/(r/m — 2)w -~t/2w as m — s- oo. 

The matrix of the second derivatives necessary to in- 
clude the Gaussian fluctuations around the saddles, has 
two distinct positive eigenvalues, r/2 and t, which are, 
respectively, r(n — r)-fold and (n — r)(n — r—l) /2-fold de- 
generate. There exist a null eigenvalue, r(r— 3)/2-fold de- 
generate, and one negative eigenvalue — r, r-fold degen- 
erate. This negative eigenvalue is responsible for the fac- 
tor (— t)~2 which generates the imaginary part for odd 
values of r. The total number of solutions that can be ob- 
tained by switching the sign of qafs is [S'(r/m)]™ 2'"/^ 
as TO ^ cxD. (This limit exists if we assume r/m is even.) 
From (|49|l . we can see that the saddles with smallest r 
dominate for pure imaginary w. Therefore, the leading 
contribution to the imaginary part of (|50|l should come 
from the first i.e. r = 1 term as we can see no reason 
why it should be excluded in this kind of replica symme- 
try breaking scheme. Thus 



±— exp 



^3 1 



48^2 

l)(n-2)4 



25 



i(n-l) 



[1 + 0{W^)], (50) 



where the upper and lower signs correspond to arg(u;) = 
— 7r/2 and arg(u') = 7r/2, respectively. 

We take the n ^ limit of ifSUI) using the fact "Cr = 
n{—lY~^ /r + 0{n'^), which can be derived from l/r(n — 
r + 1) = r(r — n)sin(7r(r — ri))/7r. We finally obtain 
the discontinuity of the imaginary part of the partition 
function ^ in the limit n ^ across the cut on the 
imaginary w (w^ < Q) axis as 



lim Im 



Z(")(ii;2) 



± exp 



48^2 



(51) 



We cannot obtain the precise prefactor to the exponential 
term because we have neglected the contributions of the 
massless eigenvalue. In principle, this soft mode could be 
integrated out by identifying the underlying symmetries 
associated with the saddle point. This does not seem ob- 
vious to us. However, its contribution is subdominant to 
that from the exponential and we shall ignore its contri- 
bution. 

The discontinuity in the imaginary part of the partition 
function can be used to extract the coefficients of the 
perturbation expansion. To do that we write the above 
result in terms of the expansion parameter 172. Recalling 
T = 2t, we can write the right hand side of Eq. H51() 
as ±exp(l/6(72) for arg(52) = zpjr. This is the same as 
Eq. H24|) up to the undetermined prefactor. Therefore we 
obtain exactly the same high-order behavior as in H25|) 



(52) 



Without the spherical model mapping one would have 
had reservations about the likely correctness of the 
replica procedure used. 



VI. HIGH ORDER TERMS OF THE e 
EXPANSION 

The starting point for obtaining the large order form of 
the e expansion is the Hamiltonian of Eq. but without 
the quartic terms: 



n 



w 



a, 13 



(53) 



Our treatment closely follows that of McKane 7] . The 
saddle points which are the analogue of Eq. H43|l of the 
toy replica calculation are the instantons which satisfy 
the equation 



= -rqap qaiQ'yp- 

7 



(54) 



For evaluating the high-order coefficients in the e expan- 
sion we can set t — and look for a solution of the form 
qai3 = w~^daf34>cir) in d = 6 dimensions. Such a so- 
lution, which decouples replica indices from the spatial 
dependence r, exists if 



da(3 — E da^d. 



7/3 > 



(55) 



and 



V20,(r) = </>2(r). (56) 
There are spherically symmetric solutions of Eq. H56|) : 

24A2 



Mr) 



"[A2r2 + 1]2' 



(57) 
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where the parameter A reflects the dilatation invariance 
of Eq. (|56|l . We shaU take for da(3 the rephca symme- 
try broken solution of scheme II. Then the energy of the 
instanton is 

a, (3 



Using the result that 



<fv (V0,(r))' = - / d^r03(^) = 



11527r-^ 



the energy of the instanton is 



E = 



11527r-^ 



48^2 



(58) 



(59) 



(60) 



where gj^ — Kqw and Kq = 5*6/ (27r)'' and 5*6 = it is the 
surface area of a six dimensional sphere of unit radius. 

The leading terms in the large order behavior of the e 
expansion are obtained by replacin g q\ by its fixed point 
value, which to lowest order is e/2 |l2l | and by the usual 
saddle-point arguments the coefficient of €^ for large K 
for any critical exponent goes like 



K 



(61) 



The next most dominant term is a factor of the form 
. The value of b depends on the critical exponent 
being studied and is beyond the scope of this paper. To 
determine its value a treatment is needed of the massless 
modes which arise in the Gaussian fluctuations around 
the instanton solution. 

Inspection of the first three terms in the e expansion 
for rj and show that these terms are not 

growing anything like as rapidly as predicted at large K. 
It is perhaps not surprising, therefore, that a Pade-Borel 
analysis of the series does not yield good numerical values 
for the critical exponents in three dimensions. 



field theory of spin glasses. By mapping this to the 
problem of critical finite-size corrections in a modified 
spherical spin glasses, we have determined the high-order 
behavior of the perturbation expansion coefficients. To 
the leading order, the coefficients alternate in sign, but 
there is a subleading contribution where the terms in 
the perturbation series show a cosine-like oscillation. In 
practice, the effects of these sub-dominant terms will be 
small, making a simple Pade-Borel resummation of the 
series useful, as was found to be the case in a similar 
situation for the disordered ferromagnet 0| . 

Non-perturbative terms are also present in spin glasses. 
These are Griffiths singularities and arise from regions 
where the values of the couplings Jy produce a smaller 
amount of frustration and hence a locally enhanced tran- 
sition temperature. A discussion of their form has been 
given in Ref. [23|. Similar singularities exist for disor- 
dered ferromagnets and it is widely believed that their 
effects are very small. To our knowledge no quantitative 
discussion of these singularities has been made for spin 
glasses and their study remains to be done. (The toy 
problem, because it is zero-dimensional, is free of Grif- 
fiths singularities). 

The e expansion for the critical exponents gives dis- 
appointing results as regards applications to real spin 
glasses. This is not just due to the fact that e = 3 in three 
dimensions as in Ref. good results were obtained for 
the exponents of the percolation problem in three dimen- 
sions from an e expansion with the same number of terms. 
We do not understand the origin of this problem. 

However, to our mind, the most significant remaining 
problem is what motivated this entire study. Namely, 
does perturbation theory (i.e. the loop expansion) work 
well in the spin glass phase or does the existence of 
"droplets" in finite dimensional spin glasses indicate that 
it fails completely? Our work does indicate though that 
perturbation theory is useful in the high-temperature 
phase. 
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